DAMPED LYMAN ALPHA SYSTEMS AND GALAXY FORMATION 
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ABSTRACT 

We examine the constraints on theories of galaxy formation that are obtained 
from observations of damped Lya (DL) systems, assuming they are gaseous 



protodisks in dark matter halos. Using the Press- Schechter formalism, we find 
that the mixed dark matter model, with Ohdm = 0.3, Oqdm = 0.65, Ob = 0.05, 
and h — 0.5, is ruled out because the number of galactic halos at z ~ 3 is too small 
to account for the total gaseous mass in DL systems, even under the assumption 
that all the gas in collapsed halos has settled into disks of neutral gas. The standard 
CDM model can account for the gas in DL systems if the bias is b ^ 2; the same 
is true for the CDM model with a cosmological constant, if b ^ 1.5 for A = 0.8. 
However, one still needs to assume that a fraction ^ 0.4 of the baryons in collapsed 
halos at z ~ 3 is in the form of neutral gas in disks. We also calculate the column 
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density distribution /(AT m ) of the DL systems, in terms of the surface density 
profiles of disks and the distribution of their central column densities. It is shown 
that the form of /(AT HI ) at the high end of column density is a diagnostic for the 
nature of DL systems. 

Subject headings: galaxies: formation - QSO: absorption lines - cosmology 
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1. INTRODUCTION 



Damped Lyct absorption systems are, at present, the best observational probe 
we have to study galaxy formation and evolution at high redshift (see Wolfe 1993). 
The observed systems are selected only because they lie on the line-of-sight to an 
unrelated quasar, rather than from any special property such as a high optical or 
radio luminosity. Thus, damped systems are probably the progenitors of typical 
galaxies. In this paper, we show that the observations of DL systems can be used 
as a constraint for models of galaxy formation. In § 2 we discuss how the column 
density distribution of DL systems can be used to infer the typical central column 
densities of the objects producing them. In § 3 we investigate the number of DL 
systems and the total gaseous mass they contain, as predicted by currently favored 
models of galaxy formation. The results are discussed in § 4. 



2. THE COLUMN DENSITY DISTRIBUTION FOR DISKS 

In this section, we calculate the column density distribution of DL systems that 
results from thin disks of HI gas. We write the surface HI density profile of a disk 

as 

E(r) = S F(r/r ), (1) 

where r is the radius on the disk, E is the central surface density, and tq is a typical 
scalelength; the function F(r/r ) is a decreasing function of r, and is normalized to 
F(0) = 1. If the normal to a disk is inclined to an angle 9 with respect to the line 
of sight, the cross section for intersecting the disk at radius r is 2irr dr cos(6>), and 
the column density observed is -/V H i = N F(r/r )/ cos(6>). For randomly oriented 
disks, cos(6>) is uniformly distributed, and the total cross section for intersecting a 



disk with column density larger than iVm is 

a(N m ) = 2n rdr cos 6>d cos 9, (2) 
Jo Jo 

where a = min [1, N F{r /ro)/AT HI ]. Differentiating with respect to iVm, we obtain: 
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where x\ is given by F(xi) = min[l, AT HI /AT ]. 

If the number density of systems of central column density between Nq and 
N + dN is n(N , z) dN 0: then the number of absorption lines found per unit 
redshift and per unit column density is 

f(N m , z) = fj ^< N ^ ( 4 ) 

where l(z) is the comoving distance at redshift z (see Sargent et al. 1980). If the 
systems have a distribution of scalelengths, then one simply needs to substitute 
in equations (3) by its average for all systems of a fixed central column density. 
When the column density Nm is higher than Nq for all the disks, then the HI 
column density distribution will be 

(see Barcons & Fabian 1987) independently of the forms of F and n. For lower 
column densities, the N m distribution depends on the profiles F(r/r ) and the 
distribution of Nq. For exponential disks with F(r/ro) = exp(— r/ro), and Nq = 
const., the HI column density distribution for Nm < N is 
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In Figure 1, we show the observed HI column density distribution, /(iVm) taken 
from Lanzetta et al. (1991) and Tytler (1987). The curves are the distribution 
function derived for exponential disks from equations (3), assuming that Nq is the 
same for all disks, for the cases N Q = 10 21 , 10 21 5 , 10 22 cm~ 2 . The amplitudes of 
the curves are fitted to the observed points. 

At present, the central surface brightness of normal galaxies are remarkably 
constant, with the value ~ 140L©pc -2 (Freeman 1970), or a column density ~ 
10 23 proton cm -2 assuming a mass-to-light ratio T ~ 5T for disks. If the DL 
systems have similar surface densities in gas, no steepening of the distribution of 
Nui should have been observed yet. But if most DL systems had much lower central 
column densities, a turnover of /(iVm) could be present in the observed range. In 
fact, if the gaseous disks at high redshift had profiles similar to the present HI disks, 
a turnover should already have been observed (see Rao & Briggs 1993). The point 
in Figure 1 at the highest column density may be an indication of such a turnover, 
and therefore of a large number of disks with No ^ 10 215 at high redshift. The 
statistics, however, are not yet conclusive. 

CONSTRAINING MODELS BY THE OBSERVATION OF DL SYSTEMS 

According to most current models, the formation of structure in the universe 
occurred through the growth of small inhomogeneities via gravitational instability. 
In this scenario, galaxies are assumed to form by cooling and condensation of gas 
in the potential wells of dark halos (e.g., White & Rees 1978). In this section 
we use the Press- Schechter formalism (Press & Schechter 1974) to estimate the 
abundance of dark halos in different models. We then combine this with some 
plausible assumptions on galaxy formation to predict the observed number of DL 
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systems and the gaseous mass they contain. 

We consider three different models for the formation of structure. Model 1 is 
the standard CDM model, with Ocdm = 0.95, Ob = 0.05 and h = 0.5 (here and in 
the following we adopt a baryon density parameter Ob = 0.0125 h~ 2 from primordial 
nucleosynthesis, where h is the Hubble constant H divided by 100kms _1 Mpc _1 , 
see Walker et al. 1991). Model 2 has a cosmological constant Oa = 0.78, Ocdm = 
0.2, Ob = 0.02 and h = 0.8; we call this model CDMA. In both of these models, we 
use the fit to the CDM power spectrum in equation (G3) of Bardeen et al. (1986). 
Model 3 (called MDM; mixed dark matter) has Ocdm = 0.65, Ohdm = 0.3, Ob = 
0.05 and h = 0.5. For this model we use the power spectrum in equation (1) of 
Klypin et al. (1993). The biasing parameters are varied for all models, and are 
chosen to bracket values currently favored by observations. These three models are 
of special interest in the "post-COBE" epoch of cosmogony (see Bond 1993 for a 
review) . 

In the Press-Schechter formalism, the comoving number density of halos as a 
function of halo mass M and redshift z can be written as 



n(M, z) dM = - — -777T 77 t _ . c — rexp 



2A(r ,z) 



where M is the mass in a top-hat window with radius r (comoving radius in 
present unit): M = 4/rpo^o/3, with po being the present cosmic mass density in 
matter. The quantity A(ro, z) is the rms mass fluctuation in a top- hat window with 
comoving radius r , linearly evolved to redshift z, and is determined by the inital 
power spectrum. The threshold 5 C is chosen to be 1.68 irrespective of models (see 
White et al. 1993 for a discussion). Assuming the spherical collapse model, one can 
define a virialized radius, r v ; r , for each halo virialized at z: 

r vir = r (l + z)-\l + 178O0 - 6 )- 1 / 3 , (8) 
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where Qo is the present-day density parameter in matter. This radius is assumed 
to define the region within which the total virialized mass of an object should be 
calculated (White et al. 1993). The circular velocity v c for a virialized halo can 
then be written as: 

1/2 

v c = (™) = 2" 1 /^ oro oJ/2 (1 + z) i/2 [x + 1780 -o. 6] 1/6 (g) 

\ ^vir / 

We shall calculate the abundance of virialized halos as a function of circular velocity 
and redshift, using equations (7, 9). 

Figure 2 shows the density of baryons (scaled to zero redshift, in unit of present 
critical density) in virialized halos as a function of redshift. This is calculated from 
equation (7) by assuming that the fraction of baryons in dark halos is the same as 
the global value in each model. The results are shown for three different intervals 
of halo circular velocities, as shown in the figures. If DL systems are due to neutral 
gas in systems which form normal spiral galaxies, then the relevant interval may be 
v c = 100 — 250kms _1 . Gaseous disks might also form in systems of lower circular 
velocity, but probably not below v c = 50kms _1 , since smaller objects would have 
a long cooling time due to the ionizing background (Efstathiou 1992). The upper 
limit corresponds to the maximum observed circular velocities in spiral galaxies, 
~ 250 km s -1 , but this is also uncertain since gaseous disks might survive for some 
time after their host halo merges with a more massive one. Clearly, the curves 
should be interpreted only as upper limits to the mass that may be present in DL 
systems, since some of the gas in halos may be ionized, or form stars rapidly, or be 
expelled from the halos in winds from supernova explosions. The data points shown 
by solid circles are the observed cosmic mass density of HI gas in DL systems, Qd, 
adopted from Lanzetta et al. (1993) for the cosmological models in consideration. 
The figure shows that the baryonic mass in virialized halos at high redshifts depends 
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strongly on the bias parameter. The large excess of the predicted at low redshifts 
is not a problem for the models, since most of the gas has probably been used up 
to form stars at present. In fact, observations of present galaxies suggest that the 
mass of neutral gas in galaxies is only ~ 10% of the total mass in stars. 

At high redshifts, the baryonic mass in galactic halos may still be larger than 
that observed in DL systems in both the CDM and CDMA models, with the bias 
parameter given in the figures. For the MDM model, however, the baryonic mass in 
virialized halos is not enough to explain the DL systems. The deficit is very large 
for the model with b = 1.5, a case favored by a normalization according to COBE 
results (see e.g., Klypin et al. 1993; Jing et al. 1993). So if DL systems are indeed 
due to gas in galactic-sized virialized halos, then the MDM model based on the 
COBE normalization can be ruled out. 

Similar conclusions have been reached recently by Subramanian & 
Padmanabhan (1994), although their values for Qd predicted by the models are 
lower than ours, because they assume a larger mass for the halos giving rise to 
damped systems. 

Figure 3 shows the comparison of the predicted rate of incidence with that 
observed. The data points are adopted from Lanzetta et al. (1993). The number of 
absorption systems above a given column density depends not only on the mass of 
gas in each halo, but also on its spatial distribution, and is consequently much more 
uncertain and model dependent than the total mass in damped systems. Here, 
we will assume that the gaseous disks causing the DL systems have exponential 
profiles with scalelengths that are the same as those of stellar disks in spiral 
galaxies at the present time, for a fixed circular velocity. Assuming also the 
Tully-Fisher relationship, and a constant central surface brightness for spiral 
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disks, the scalelength is Ro = 3.5(f c /220kms~ 1 ) 2 /i~ 1 kpc; the normalization is 
chosen to give the right range of scalelengths observed (van der Kruit 1987), for 
100 < v c < 250 km s -1 . We then use equation (4), integrated over iV HI , to find 
the total number of systems per unit redshift above a given column density. We 
take A^hi = 10 20 3 cm -2 for the minimum column density required to identify an 
absorption line as a damped system (see Wolfe 1993). 

From the results shown in Figure 3, we see that the predictions of the CDM 
model with b = 1, 2, and the CDMA model with 6=1 may be consistent with 
observations, if other effects do not reduce the cross section significantly. In the 
CDMA model, n(z) at high z depends strongly on the bias parameter b. If b « 
1.5 (see Cen, Gnedin & Ostriker, 1993), the CDMA model may have difficulty in 
explaining the observed n(z). But for the MDM model with b > 1, there is a much 
clearer disagreement with observations. The discrepancy at high redshift in this 
model is larger in n(z) than in Qd(z) (see Fig. 2), because the average HI column 
density predicted by the model is larger than observed. This is primarily due to the 
fact that the central density of disks goes as v~ x in our disc model and that most 
mass is in small haloes at high redshifts (see Fig. 2). 



4. DISCUSSION 



In this paper, we have seen that the large amount of mass in neutral gas 
observed in DL systems at high redshifts rules out the MDM model. The reason 
is that the fraction of gas in collapsed halos which can produce DL systems is too 
small. Even in the CDM and CDMA models, one needs to assume that a large 
fraction of the gaseous mass in galactic halos at z ~ 3 is in disks of neutral gas. 
Our calculated number of DL systems is more uncertain, since it depends on the 



radial distribution of the gas in the disks. The number of DL systems could be larger 
if the radius where the column density is equal to 10 20,3 cm -2 (the threshold for 
DL systems) was increased. This could be done while keeping most of the gaseous 
mass close to the center, where it would produce higher column density systems, 
and would turn into the stellar disks in spiral galaxies at present. The lower column 
density regions of the disks at large radius might be destroyed before forming many 
stars, due to supernova explosions, or infalling material resulting in shock-heating 
of the gas. In fact, most HI disks in present galaxies are observed to flare at large 
radius (van Gorkom 1993 and references therein). 

A second possibility is that there is a population of low column density disks, 
which extend to large radius. The star formation rate in such disks might be much 
lower than in regular spiral galaxies, and they might either be destroyed when the 
galactic halos where they reside merge, or might evolve into low surface brightness 
galaxies at present, which could be difficult to detect (McGaugh 1994). 

Further constraints on the nature of DL systems may come through a 
better determination of the column density distribution (which depends on the 
surface density profiles and the maximum column density of the disks), and from 
measurements of their transverse size, either in gravitationally lensed quasars, or 
by direct imaging of the absorbing galaxies (Briggs et al. 1989; Wolfe et al. 1992). 

We would like to acknowledge Stacy McGaugh, T. Padmanabhan, Martin Rees 
and Simon White for helpful discussions. HJM acknowledges support from a SERC 
Postdoctoral Fellowship. JM thanks SERC for support in Cambridge, and the W. 
M. Keck Foundation for support in Princeton. 
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FIGURE CAPTIONS 



Fig. 1: The HI column density distribution function for damped Lyman alpha 
systems. The five observational data points at high A^hi are adopted from Lanzetta 
et al. (1991); the other one is adopted from Tytler (1987). The curves (discussed 
in subsection 3.1) show the predictions of a model in which all damped Lyman 
alpha absorbers are exponential disks with central column densities iVo shown in 
the figure. The amplitudes of the curves are adjusted to the observation. 

Fig. 2: The density of baryons (scaled to zero redshift,in unit of present critical 
density) in virialized halos as a function of redshift, predicted by CDM (left panel), 
CDMA (middle) and MDM (right). For each model, results for two bias parameters 
are shown. For each bias parameter, results are shown for three different intervals 
of halo circular velocities: v c > 50, 50 — 250, 100 — 250 km s -1 , with the higher curve 
corresponding to the case with a larger interval. The data points for the DL systems 
(solid circles) are adopted from Lanzetta et al. (1993), adjusted according to the 
cosmological model in consideration. The upper cross shows the density parameter 
of baryons in stars of present-day galaxies. The lower cross shows the HI mass in 
current spirals given by Rao & Briggs (1993). 

Fig. 3: The rate of incidence, n(z), of damped Lyman alpha absorptions 
predicted by CDM (left panel), CDMA (middle) and MDM (right). The gas mass 
in halos is assumed to settle in exponential disks with scalelengths given by the 
Tully-Fisher relation (see text). 
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